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تحلیل پاسخ دینامیکی سازه های بلند تحت تأثیر نیروی محوری * 


() 0) 


مهرداد محمدنژاد حسن حاجی کاظمی 
چکیده _ در این تحفیق, معادلات انتگرالی تضعیف شده جهت محاسبه ف رکانس های طبیعی ارتعاش سازه های بلند دارای جرم و سبحتی 
متغیر در ارتفاع و تحت تأثیر نیروی محوری ارائه شده است. معادله دیفرانسیل حاکم بر ارتعا شآزاد یک تی رکنسول غیرمنشوری با سختی 
برشی و سختی خمشی و تحت تأثیر نیروی محوری, از طریق چهار مرتبه انتگرال گیری پی در پی به فرم تضعیف شده ‏ یآن تبدیل شده 
است. با تقریب تابع مد شکل ارتعاش توسط یک سری توانی, معادله انتگرالی به یک دستگاه معادلات جبری عطی تبدیل شده است. در 
نهایت با محاسبه جواب غیر بدیهی دستگاه معادلات » فرکانس های طبیعی ارتعاش سازه بلند بدس ت آمده است. سازه مورد نظر در نرم افزار 
2-0 مدل سازی شده و نتایج تحلیل به کمک روش ارائه شده, نتایج مدل سازی در نرم افزار و همچنین نتایج بدس ت آمده از 


روشهای سایر محففین با یکدیگر مفایسه شده است. 


واژه های کلیدی سازه بلند. تیر برشی- خمشی, فرم تضعیف شده معادله انتگرالی. نیروی محوری» فرکانس طبیعی ارتعاش. 
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*تاریخ دریافت مقاله ۹۵/۳/۱۸ و تاریخ پذیرش آن ۹0/۵/۱ می‌باشد. 
(۱) نویسندهٌ مسئول. دانشجوی دکتری سازه. گروه مهندسی عمران دانشگاه فردوسی مشهد. را ۱۰۱۱۰۱۱۳۵۱۸ 


(۲) عضو هیات علمی گروه مهندسی عمران دانشگاه فردوسی مشهد. 
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مقد مه 
رفتار تعداد زیادی از پدیده‌های طبیعی را می‌توان به 
کمک معادلات دیفرانسیل معمولی یا مشتق جزئی 
مدل‌سازی کرد. به عنوان مثال ارتعاش تیرهای الاستیک 
با خصوصیات جرم و سختی متغیّر به کمک معادلات 
دیفرانسیل مشتق جزئی با ضرایب متغیر مدل می شود. 
روشهای تحلیل سازه بلند را می توان به دو دسته کلی 
: روشهای تقریبی و ۲: روشهای تحلیل اجزا محدود 
دسته بندی نمود. جهت تحلیل تقریبی یک سازه بلند. از 
یک تیر کنسول با مقطع متغیر به عنوان جایگزین سازه 
اصلی استفاده می شود. مدلسازی سازه بلند توسط یک 
تیر جایگزین به کمک تئوری تقریبی تیر جایگزین 
(صههها اممجصمع۵۵۱2 ) انجام می شود. با حل معادله 
حاکم بر ارتعاش این تیر جایگزین فرکانس های طبیعی 
ارتعاش سازه اصلی و توابع مد شکل محاسبه می شوند. 
بنابراین تحلیل تير جایگزین به جای سازه اصلی دارای 
ساده سازی ها و تقریب هایی می باشد. روشهای تحلیل 
اخا موه بر باه لادم کته انار اه این 
روش ها نیاز به حل همزمان هزاران معادله خطی به 
منظور ارائه نتایج عددی و کاربردی تحلیل سازه دارد. 
همچنین مدلسازی یک سازه بلند در نرم افزار تحلیل 
اجزا محدود بدلیل داشتن المان های سازه ای زیاد بسیار 
یف تا تن و ی اس ری 
برای مراحل نهایی و تحلیل دقیق یک سازه بلند مناسب 
زا اتف تال ایتک کهروتمای تقرس ای 
سازه, به دلیل سادگی و سرعت عمل بالاء در مراحل 
وله طراشی .و فان مدای شکارم هار 
روش ارائه شده در این تحقیق یک دید روشن از تحلیل 
ارتعاش آزاد سازه بلند ارائه می کند. سادگی و دقت 
بالای روش ارائه شده در این تحقیق موجب می شود 
کی لین تفریی وراج معدمانی طراجی ,یک 
سازه بلند ( که نیازی به مدلسازی دقیق نرم افزاری 
وجود ندارد) قابل استفاده باشد. علاوه بر این: جهت 
بررسی میزان تأثیر پارامترهای مختلف سازه ای بر 
مشخصات استاتیکی و دینامیکی یک سازه بلند. استفاده 
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تحلیل پاسخ دینامیکی سازه های بلند تحت تثیر نیروی محوری 


از روشهای تحلیلی راحتتر و سرراست تر می باشد. با 
توجه به زمان اندک مورد نیاز جهت مدلسازی و 
محاسبات در روش ارائه شده ‏ این روش در تحلیل 
تقریبی سازه مفید واقع می شود. نتایج تحلیل تقریبی 
می تواند در جهت کاهش خطاهای مدلسازی دقیق 
سازه در نرم افزار تحلیل اجزا محدود. نخمین صحیح 
ابعاد اولیه المان های سازه ای و مشخصات مصالح 
مورد بررسی و ... مفید واقع گردد. 

تاکنون تحقیقات زیادی در خصوص مدلسازی و 
تحلیل تقریبی سازه های بلند ارائه شده اتبتا: در سال 
۳۰۱۲ قزاق و همکاران [1] یک توری اصلاح شده 
برای تفکیک یک سازه بلند دارای سختی خمشی و 
برشی به دو زير سازه مجزا (یکی برشی و دیگری 
خمشی) ارائه کرده اند. در سال ۲۰۱۶ پارک و همکاران 
اک ان ام تا گم شا نها 
بلند دارای سیستم سازه ای قاب-دیوار برشی ارائه کرده 
اند. یک سیستم قاب-دیوار از یک دیوار برشی با رفتار 
خمشی به عنوان هسته و یک قاب خمشی با رفتار 
برشی تشکیل می شود. در سال ۳۰ هکل و 
هت‌کاران زارف کانیی‌های: یی یک نموه شک 
از قاب لوله ای و دیوار برشی را به کمک روش 
گالرکین و توابع عذا3-5 محاسبه کرده اند. در این 
تحقیق, معادله حاکم بر ارتعاش آزاد یک سازه بلند با 
به فرم ماتریسی آن تبدیل شده است. در سال ۲۰۱۲ 
کامگار و همکاران [4] فرکانس طبیعی یک سازه بلند 
مرکب از قاب لوله ای هسته برشی. کمربند خرپایی و 
مهاربازویی که دارای ناپپونسکین هایی در مقطع عرضی 
می باشد را محاسبه کرده اند . در این تحقیق فرض شده 
است سختی و جرم سازه به صورت پلکانی و ناگهانی 
در ارتفاع سازه کاهش می یابند. از این رو سازه بلند به 
8 قسمت در ارتفاع تفکیک شده و برای هر قسمت 
سختی و جرم ثابت فرض شده است. با اعمال شرایط 
مرزی در محل فصل مشترک هر دو قطعه پی در پی 1 
اه هبدن نا راکب تایت ار ضاظر نا .9 شمیت 


سال سی و یکم شمارة دوه ۱۳۹۷ 


مهرداد محمدنژادسحسن حاج ی کاظمی 


تکیک: که تعاضا شله اشتته با حل هریاد لته 
روش سری توانی (5010000 56165 00۷6۲ ) فرکانس 
قع طاقن اوه تا سایق سنا بر 
سالهای ۲۰۱۲ و ۲۰۱۶ ملکی نژاد و همکاران [5-6] 
روشهایی بر مبنای روابط انرژی برای محاسبه فرکانس 
های طبیعی و توابع مد شکل ارتعاش سازه های بلند 
دارای سیستم قاب لوله ای. هسته برشی و کمربند 
خرپایی ارائه کرده اند. در این تحفیق. سازه قاب 
محیطی به کمک یک تیر کنسول با مقطع تیر-جعبه و 
سیستم مرکب کمربند خرپایی-مهار بازویی به کمک 
یک فنر دورانی مستقر در مکان کمربند مدلسازی شده 
انسکه هن سال ۲۰۹۱۳ الک زر روانظ بزای ی 
پیچشی سازه های بلند نامتقارن در پلان با سیستم باربر 
جانبی قاب خمشی, دیوار برشی و هسته مرکزی ارائه 
کرده است. در سال ۲۰۱۵ محمدنژاد [8] و در سال 
۶ محمدنژاد و همکاران [9] فرکانس های طبیعی 
ارتعاش تیرهایی با خصوصیات جرم و سختی متغیر در 
طول تیر و تحت تأثیر نیروی محوری و اینرسی دورانی 
را به کمک معادلات انتگرالی تضعیف شده محاسبه 
کرده اند. در اين تحقیق, معادله حاکم بر ارتعاش 
تیرهای مختلف از طریق انتگرال گیری های پی در پی 
به فرم تضعیف شده آن تبدیل شده است. در سال ۲۰۱۵ 
صفاری و همکاران [10] فرکانس ارتعاشی سازه های 
بلند دارای سختی خمشی و برشی را به کمک معادلات 
انتگرالی تضعیف شده محاسبه نموده اند. در این تحقیق, 
سازه بلند با رفتار برشی و خمشی به دو سازه مجزا 
یکی با رفتار حمشی و دیگری با رفتار برشی تفکیک 
شده است. معادلات حاکم بر ارتعاش تیر برشی و تیر 
خمشی به کمک معادلات انتگرالی تضعیف شده حل 
شده و فرکانس های ارتعاش برشی و ارتعاش خحمشی 
ابا تایه تم اس فر خاش ای انیب فیک 
ترکیب فرکانس های برشی و خمشی محاسبه شده 
اه رای ۱۳۱۲ مار و هار ان ۱ کش 


های طبیعی ارتعاش و بار کمانشی تیر برنولی با 


سال سی و یکم شماره دی ۱۳۹۷ 


3 


خصوصیات سختی و جرم متغیر در طول تیر و تحت 
اثر نیروی محوری متغیر را به کمک معادلات انتگرالی 
تضعیف شده محاسبه کرده اند. در اين تحقیق فرکانس 
های طبیعی ارتعاش و بار کمانشی تیرهایی با شرایط 
تکیه گاهی مختلف محاسبه شده است. تحقیق های 
زیادی در زمینه تحلیل ارتعاش آزاد سازه های بلند به 
کتک رهشیان تلف ار ایه نات رف ایادخ 
تعداد کمی از این تحقیقات ار نیروی محوری بر روی 
رفتار ارتعاشی سازه بلند بررسی شده است. از آنجایی 
که نیروی محوری موجب کاهش سختی سازه می گردد 
در نظر گرفتن این پارامتر در محاسبه فرکانس طبیعی 
ارتعاش امری ضروری به نظر می رسد. 

در تحقیق پیش رو فرکانسهای طبیعی ارتعاش 
یک سازه بلند دارای جرم و سختی متغیر در ارتفاع 
فا نان کتروع: مزر مخاسمه شم ات عادنة 
حاکم بر ارتعاش از طریق چهار مرتبه انتگرال گیری پی 
در پی به فرم تضعیف شده آن تبدیل می شود. اکثر 
روشهای تحلیلی ارائه شده تاکنون جهت تحلیل ارتعاش 
آزاد سازه بلند. سختی و جرم سازه را در ارتفاع 
یکنواخت فرض کرده اند. در روش ارائه شده در این 
تحفیق. خصوصیات سختی و جرم سازه در ارتفاع 
متغیر فرض شده است. که این نوع مدلسازی به واقعیت 
موجود در سازه ها نزدیک تر می باشد. همچنین. تحلیل 
ارتعاش آزاد یک تیر با مقطع ثابت (که متناظر با مدل 
یک سازه بلند با خصوصیات سختی و جرم ثابت می 
باشد) به مراتب ساده تر از تحلیل یک تیر با مقطع غیر 
منشوری (که متناظر با مدل یک سازه بلند با 
خصوصیات سختی و جرم متغیر است) می باشد. در 
اکثر روشهای ارائه شده تاکنون سازه بلند که دارای یک 
محیط پیوسته با بینهایت درجه آزادی می باشد به لا 
سازه یک درجه آزادی معادل, متناظر با آا مد ارتعاشی 
تبدیل شده و برای هر مد ارتعاشی یکبار معادله حاکم 


حل شده و سپس نتایج تحلیل بایکدیگر ترکیب شده 
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اند. دو اپراد اساسی به این ایده وارد است. نخست 
که اس بارس اک مخ یشرت وین هن 
است و دوم اينکه نیاز به لا بار حل معادله حاکم بر 
ارتعاش می باشد. در حالی که در روش پیش رو به 
کمک مدلسازی پیوسته سازه بلند ( مصسصتاجم 
00108), که به واقعیت سازه نزدیک تر است. تنها 
یک‌باز معادله سخاکم بر ازتفاشن سیم پپوسته نحل شاه 
ای منم هر انظ( ۱۵ تاک سا 
مد ارتعاشی که مد نظر باشد دسترسی داریم و نیازی به 
ترکیب جواب ها نیز نخواهيم داشت. تاکنون جهعت 
تبدیل معادله حاکم بر ارتعاش آزاد سازه بلند به یک 
معادله قابل حل. روشهای حل بسیار طولانی همراه با 
بات پاش هدزاش ات هرانک ور 
این تحقیق تنها از طریق ۶ مرتبه انتگرال گیری پی در 
پی از معادله حاکم به یک معادله قابل حل دسترسی پیدا 
کرده انش که ای ما ناد کی زرشن یی وتان 
می دهد. در معدود روشهای ارائه شده تاکنون برای 
تحلیل ارتعاش آزاد سازه هایی با خصوصیات سختی و 
جرم متغیر در ارتفاع» روش حل بدین صورت می باشد 
که سازه غیر منشوری به 9 قسمت تقسیم شده و فرضص 
شده است سختی و جرم سازه در هر قسمت ثابت می 
باشد و شرایط مرزی بین هر دو قطعه پی در پی ارضا 
ید ان فا یزاین فاوبه سل دنه دیق انسیا سا 
ضرائب ثابت متناظر با 9 قسمت سازه وجود دارد. در 
حالی که در روش تحقیق پیش رو نیازی به این تقسیم 
بندی وجود نداشته و برای تغییرات سختی و جرم سازه 
در ارتفاع محدودیتی وجود ندارد و تنها نیاز به بررسی 
شرایط مرزی در پای سازه و بام سازه می باشد. که اینن 
اسان سا شتا له وا ای کش دای 


محاسبه فرکانس طبیعی ارتعاش سازه بلند 


سازه بلند توسط یک تیر کنسول باسختی خمشی. 


نشرية مهندسی عمران فردوسی 


تحلیل پاسخ دینامیکی سازه های بلند تحت تثیر نیروی محوری 


سختی برشی و جرم واحد طول متغیر در ارتفاع تحت 
تأثیر نیروی محوری ناشی از وزن سازه مدلسازی می 
شود. توابع سختی و جرم تیر جایگزین به کمک روابط 
تقریبی محاسبه می شود. معادله حاکم بر ارتعاش یک 
تیر برشی -خمشی(0620 9۳6۵۲-۳۴600108) غیر منشوری 
تحت تأثیر نیروی محوری در نظر گرفته شده و به 
کمک معادلات انتگرالی تضعیف شده فرکانس های 
طبیعی ارتعاش این تیرکه در حقیقت فرکانس های سازه 
اصلی می باشند. محاسبه می شود. تیر برشیسخمشی 
یک حالت خاص از تیر تیموشنکو می باشد که در 
معادله حاکم بر ارتعماش آن از اثرات اینرسی دورانی 
صرفنظر شده است. در سال ۳ کامگار و رهگذر 
[23] در تحقیقی نشان داده اند که اثر اینرسی دورانی بر 
فرکانس طبیعی سازه بلند ناچیز می باشد. 


تبدیل معادله حاکم به فرم تضعیف شده. با استفاده از 
روش انرژی و اصل همیلتون (2۳001016 و ممنانصه1) 
معادله حاکم بر ارتعاش سازه بلند با در نظر گرفتن 
اثرات تغییر شکل برشی, تغییر شکل خمشی و نیروی 


محوری به صورت زير محاسبه شده است [4: 


]۱00۷۵ | 0 0۷ بل 
ی 


0)26,1( 2 0, >> 
(۱) 

در رابطه (۰6۱ (۷۷)26,۸ ,(26) ویک ب(26) وکا ,(10)26 
,(11,۱۷)۸6 و (406,1 به ترتیب تابع جابجایی جانبی 
سازه تابع سختی خمشی. تابع سختی برشی, تابع جرم 
واحد طول سازه تابع نیروی محوری سازه. ارتفاع سازه 
و تابع بارگذاری جانبی سازه می باشند (شکل ۱). 
جهت بررسی ارتعاش آزاد 20 (9)2,1 در نظر گرفته 


می شود. 


سال سی و یکم شمارة دوه ۱۳۹۷ 


مهرداد محمدنژادسحسن حاج ی کاظمی 


۳۳ 
7 4 
۳ 
۱۱ )و1 
سچ ۳ 
۳ 
۳ 
۳ 
)۷ 
(الف) 0 (ب) 


شکل ۱ مدلسازی تقریبی سازه بلند توسط تیر جایگزین تحت اثر 
نیروی محوری. الف) سازه اصلی. ب تیر جایگزین سازه اصلی 


با فرض یک ارتعاش هارمونیک برای سازه. تابع 
جابجایی جانبی تير به صورت زیر در نظر گرفته می 
شود: 
9 و() ۷ < (2,۱) ۱۷ 
در رابطهٌ (۲ (۷)2 و 2 به ترتیب تابع مد 
فک ناشن و فرکاشی ظیی. اتقعاشن ,سازورای 
معادله تک متغیره بر حسب پارامتر مکان به صورت 


زیر بداست می آ بل 


و[ ون 
,0-(06) ۷ (6) م تج سر ك )2 ۳ ك 
>> 


(۳ 


به منظور سهولت محاسبات. تغیر متغیر 3 » در 

نظر گرفته شده و در رابطه (۳) جایگذاری می شود: 
0 رم 6 1 2 اک 
۳ ۱ ان ب 0 
ر[_ 3 )۳( 1 0 ۳3 )و 23 
9 1 
0 < () )0 62*۲۲ + ها با 
۹ ۹ 
1> > 0 


(1) 


معادله (۶) معادله ازتعاش آزاد سازه بل بر خسب 


سال سی و یکم شمارة دو ۱۳۹۷ 


2 


به فرم تضعیف شده آن. از طرفین معادله (6) دو مرتبه 
یت یازا در بان 0 شا و اعکرآل گیزی :مس 


شود. نتیجه, معادلات انتگرالی زیر می باشند: 
6 2 اک 2 
() ۱۷ ()۱ ۲۲۰ + () ۱۷ () یک1 ۲۲ 
۹ 1 
1 1 
تهب 


بل - همه رت 


2 
)۱۱ ۲۲ +( ۷( :۳۲1 
ی ِِ » 
(۱۲۵ ۲ - (9)یک 11و 
)+6 < ریخ 
۹ 
0( 


همچنین دو مرتبه انتگرال گیری از معادله () 
نسبت به پارامتر « در بازه 0 تا ۶ معادلات انتگرالی زیر 


را نتیجه می دهد: 


- (۳(۷)۵ + (۵) ۳ | 
-(8) 1 (9-ک) ۲۲ . 

5 | + (و) 162 - (9) ۸ (9- 0 112 
(0)5 *(-) 1 


- 00۷00 16 + )رک - 


بل + در + رگ 


۹2 
2 
< ()۱۷ () و1 - 3(5) ۱۷ (5,) ۳۳ 
0 
بل +0 + شرت یز 1 
2 6 

)۸( 
در رابطه (۸) تابع (9,«) ,طبه صورت زیر محاسبه 

می شود: 


نشریه مهندسی عمران فردوسی 


31 
+(و) یک (-) ۷2 - (وب) رط 
۵ ارم( (و ها 
-(۵) 2 +(۵) ۲ وس[ 
هه ارو بر مک (و- یم 
(٩)‏ 


معا ۳( فرع مش کم معادله: عطاکم رشن 
ارتعاش سازه بلند با در نظر گرفتن اثرات نیروی 
محوری می باشد. در معادله (۸) پارامترهای 
رن« و),و) ,ی ثابتهای انتگرال گیری می باشند که به 
هورق سانش و شتااسی تاک ار مات 
سازه بلند محاسبه می شوند. در قسمت بعد. از معادلات 
(۵) الی (۸) جهت تعیین ثابت های انتگرال گیری 
تاو یقرت 


شرایط مرژی. شرایط مرزی اساسی در تکیه گاه تیر به 


صورت زیر در نظر گرفته می شود: 
2-0 (۷۷)0 20 ۲ 


(۷۰) و (۱۱) 

شرایط مرزی هندسی بسته به نوع سیستم باربر 
پارامترهایی نظیر نیروی محوری, اینرسی دورانی» تغییر 
شکل های برشی و 


مرزی هندسی ف تا واه نت ابر رون 


رت رم -لاشق: ز اانقط ق ابظ 


محوری به صورت زد پر در نظر گرفته می شود: 
6 2 
() ۷۷ () وک 11- 21 ۶ 
۹ 
0 2 
9( 1 


مس 9 1۳ 7 


۳ 
20 (0 ۷ 7 () وکا 1< 7 
۹ 


(۱۲) و (۱۳) 


نشرية مهندسی عمران فردوسی 


تحلیل پاسخ دینامیکی سازه های بلند تحت تثیر نیروی محوری 


رابطه (۱۲) مربوط به وضعیت نیروی برشی در 
انتهای آزاد تیر تحت تأثیر نیروی محوری می باشد. در 
این رابطه. ترم سوم همان نیروی برشی کلاسیک تیر 
می‌باشد ترم دوم اثر نیروی محوری و ترم اول اثر 
سختی برشی بر روی نیروی برشی تیر می باشند. 
همچنین رابطه (۱۳) مربوط به وضعیت ممان خمشی در 
انتهای آزاد تير می باشد. با اعمال شرایط مرزی (۱۰) 
الی (۱۳) در روابط (۵) الی (۸ ثابت های انتگرال 
گیری محاسبه شده و در رابطه (۸) جایگذاری می شود. 
نتیجه, معادله انتگرالی زیر بر حسب تابع مد شکل 
ارتعاش ()7 می‌باشد: 


1 ۹ 
- 5 (5) ۱ (5,<) ]+5 (۹(۱)۵,) ۳ 
0 0 
20 (۷۷)۶() وک 
)۱ 


در رابطه (۱۶) تابع (2)5,5 به صورت زیر 


محاسبه می شود: 
24 


21 
< (۹,) روط 


0 
7( زر - () 0 


022 
ار سب 
(5) 2 2 


(۱۵( 
در رابطه (۱۵) توابع (8),عو (8)ع و پارامترهای 
,0و ول به صورت زیر محاسبه می شوند: 


- (1>۵)8 ۲۳ - (و)صرو-) ۲۱۴ 2 ع (ع) رم 
(۳5 02711 - (5) ۲۲7 


9 
+(8) >( ی (0) ,8 
(۵) 14 *(و 0 (۱-9(۱۷۵) 
9 2 ِِ 
ن ان 0 0 
و 
9 
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(1) و21 ۳ 
() 21 +۱0 ۲+(10 ۲ ۲ 
۲۵۵ ۲ () 21 +(1,)0 ۲۱ (1) 2 
(1) ,21+( 1۲۷۰ + )1۱1 ِ 


(۱۸) 


تیر کنسول با جرم متمرکز در انتهای آزاد. در این 
قسمت فرض می‌شود که یک جرم متمرکز در انتهای 
تیر کنسول وجود دارد. این فرض جهت تحلیل 
سازه‌های بلند که دارای جرم متمرکز قابل توجه در 
نتهای آزاد خود می‌باشند (نظیر مخازن اب). قابل 
تتفاه: میباشتت قوش یل مشاه یت قز 
ش‌باشند:نا این فقاوت که فر این سالت: تبروی پرشی, ذر 
نتهای تیر کنسول و در محل جرم متمرکز برابر صفر 
نیست. سایر شرایط مرزی تغییر نمی‌کنند. با صرف نظر 


ز اثرات سختی برشی. در این حالت شرایط مرزی به 


صورت زير تعریف می‌شوند [24: 


 < 0, ۱۷ )0( < 0,‏ 
"1 
20 (0)-س- ,0 < 
)00 1 َ 
اک 2 
()۱۷-- () "۳1 - بر  <‏ 
1 


س ط 4 
(1) ۱۷1622۷- - ی مد 
1 


( جر 


۱۹( 
در روابط (۱۹) ۷ جرم متمرکز در انتهای آزاد تیر 
کل پاش با اعتان کرایط مروي فرن خات‌هاي 
انتگرال گیری معادله(۸) تعیین می‌شوند. با جایگذاری 
ثابت‌های انتگرال‌گیری در معادله(۸) یک معادله انتگرالی 

به شکل زير بدست می‌آید: 
)٩(6‏ ۱۷( ,<) ۳ + 5(05) ۱۷ (5,) ۳ 
0 0 
20 ()۱۷() وک[ 
(۳۰( 


سال سی و یکم شماره دی ۱۳۹۷ 
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در معادله (۲۰) تابع (9,ع)مط به صورت زیر 


محاسبه می‌شود: 
- 0 | دشک تب | ودرا 


2 04 
8 كگن 
)8( ِ 


)۳۱( 
در رابطه فوق: 
<(و) ویک( -1) + (1-9(()۵) ۲12 - (و) بع 


99*6 
6 


- (9) و21 - (65 *(1-9) 


7 
۳-۱-۳۱۷۸ 
+(9) 3117+ () بع6 < (9) بت 


(۳2)6 1 0 232۳ ِ )1( وله - (31۳1*(0 
۷ 


(۳) 


۹98 1 
1 
(3۳1*(۷0- 6160-2167117 
(1) و216 _ 
و ۳ 
۵۷۲+ (1) و616 ِ 
(1) و181 (911*(۲)۱ + 6۱۵2۲۲۲ 


:را 


5 
(۳۳ 


در نظر گرفته شده اند. تابع (5,ع) بط از رابطه 69 


انیب مین کرد: 


تبدیل معادله انتگرالی به دستکاه معادلات جبری 
خطی. ‏ تنها مجهول معادلات انتگرالی (۱4) و (۲۰) 
تابع مد شکل ارتعاش (۷69می باشد. این تابع توسط 
سری توانی زیر تقریب زده می شود: 


۳ ۲( 


9 
(ا 
در این سری ,) ضرایب مجهول سری هستند و ۲ 

بستگی به مقدار ۴ انتخاب شده دارد. با در نظر گرفتن 
یک ۳ بزرگتر تعداد جملات بیشتری از سری (۲۶) در 


نشریه مهندسی عمران فردوسی 
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نتایجم می شود. با جایگذاری سری (۲۶) در معادلات 


انتگرالی (۱6) و(۲۰) و ضرب طرفین معادله در ۳و 
انتگرال گیری از معادله حاصل در بازه 0 تا 1 نسبت 
به پارامتر ۰ س ام معادلات جبری خطی زیر بدست 
1 
0 < تا( 9 بحط) ر1 + (۲یحط) ,۳۱ + (1 ,0( 
۲<0 
طو.. و2 وکر0 < 10 
(۲۵( 
در رابطه (۲۵) توابع (تیص)ت (۲یه) 1 و 
(0,1) 112 به صورت زیرمحاسبه می شوند: 


1 
رک ود [ -- ,00 
0 


1 
یل و(وی) رط 0 [ - (1,) ,11 
00 


1 1 
وی رط 0 [ > (0,۲) و11 
00 


(۳۹0 


سیستم دستگاه معادلات جبری خطی (۲۵) را 
میتوان به فرم ماتریسی زیر بازنویسی نمود: 
۳۷ مریب01< مریبیلکا ریبی ورب[ ش] 
در رابطه (۲۷) [۸]ماتریس ضرایب و [)]ماتریس 
مجهولات می باشند. تنها پارامتر مجهول در ماتریس 
[ ]۰ فرکانس طبیعی ارتعاش سازه می باشد. با محاسبه 
جواب غیر بدیهی دستگاه معادلات(۲۷) که با مساوی 
صفر قرار دادن دترمینان ماتریس [۸]بدست می آید. 
معادله فرکانسی ثیر محاسبه می شود « ویشه های, معادله 
فرکانسی, فرکانس های طبیعی ارتعاش سازه اصلی می 
باق شاحله زر کاتشی یک نانمس با اش او درخ 
(2)۳+1می باشد بنابراین با در نظر گرفتن ۴بزرگتر 
تعداد ریشه های بیشتری (فرکانس های مدهای بالاتر) 
قابل محاسبه می باشد. همچنین با در نظر گرفتن ۲ 
بزرگتر دقت نتایج تحلیل افزايش می یابد. 


نشریة مهندسی عمران فردوسی 


تحلیل پاسخ دینامیکی سازه های بلند تحت تثیر نیروی محوری 


تابع مد شکل ارتعاش 
پس از محاسبه فرکانس های طبیعی ارتعاش سازه به 
کمک روش تشریح شده در قسمت قبل. می توان توابع 
مد شکل ارتعاش را محاسبه نمود. از آنجایی که تابع مد 
شکل ارتعاش به کمک سری توانی زیر تقریب زده شده 
است: 


۲ 
۱۷ )( < 2 (۳۸) 


جهت محاسبه تابع مد شکل ارتعاش نیاز به 
محاسبه ثابت های 7 مت می باشد. پس از محاسبه 
فرکانس طبیعی ارتعاش مد (ذ) ام :6۵ و جایگذاری آن 
در دستگاه معادلات جبری خطی (۲۷) دستگاه معادلات 
زو لیس هی اه 
(۳۹( 0< ررربم[ وأربمییم[ ن۸] 
در رابطه (۲۹) 21 ()) در نظر گرفته شده و 
جایگذاری می شود. نتیجه دستگاه معادلات جبری زیر 
می باشد: 
۳ ونم لها 


در رابطه (۳۰) ماتریس مای وریبم[ زظ]: 
رم ] «ررربم[ن۸] به صورت زیر محاسبه می 


:| «رببم[ :ظ] 


۳ / 


۳ / 


(۳۱ 
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با حل دستگاه معادلات (۳۰) ثابت های ) 


شا هو از ی 


نیروهای داخلی 
در آنالیز مودال سازه هاء برای هر مد ارتعاش یک بار 
معادل تعیین می شود. نیروی جانبی معادل در مد (1) ام 
به صورت زير محاسبه می شود(شکل 25[)۲]: 


1 
«00 0 

زمر 26(5) :22)2(۷۷ ۳ 
۳ 
0 

0 
1 
«00 )( 

نم 5( :۷۷ ( )22 ۲ < () :1 
() ۱00 
0 


(۳۲) 
در رابطه (۳۲) م5 شتاب طیفی متناظر با مد 1 ام 
ارتعاش (که وابسته به زمان تناوب ارتعاش, میرائی سازه 
و شتاب ماکزیموم زمین می باشد) و 11 ارتفاع سازه می 
باشد. با انتگرال گیری از رابطه (۳۲) بین 26 و 31 (يا بین 
#و ۱) نیروی برشی موثر بر مقطع عرضی سازه در مد 1 
ام بدست می آید (شکل ۲). نتیجه انتگرال گیری از 

رابطه (۳۲) به صورت زیر می باشد: 


11 1 
۷00 - [0 ۰ ۷0۵ 1] 1:۵ ۳۳ 
4 ۹ 


با قرار دادن 20 در رابطه (۲۳ نیروی برشی کل 
و با اتقو زو 
محاسبه می شود. لنگر خمشی پای سازه به صورت زیر 
محاسبه می شود (شکل ۲): 


1 1 
۷ ۳ 0۲ ۷: < ۴ | ۳۶( 
0 0 


سال سی و یک شماره دوب ۱۳۹۷ 


۷ 


شکل ۲ بار جانبی معادل ۰80 نیروی برشی موثر بر مقطع 
عرضی سازه ۷:۵ و نیروی برش پایه :۷ و لنگر خمشی پای سازه 
:۶ در مد( 1) ام ارتعاش 


مثال های عددی 
سازه بلند با سیستم دیوار برشی. در این مثال ۵ 
فرکانس ابتدایی یک سازه بلند ۲۷ طبقه با ۷۱ متر 
ارتفاع و جرم متمرکز در انتهای آزاد محاسبه می شود. 
سیستم باربر جانبی این سازه دیوار برشی با مقطع متغیر 
در ارتفاع سازه می باشد. از آنجایی که سیستم باربر 
جانبی تنها دارای سختی خمشی می باشد این سازه 
توسط یک تیر کنسول خمشی با جرم متمرکز در انتهای 
آزاد سل ی شود از انرات برش و آزسی قوراین 
صرف نظر شده است. تابع نیروی محوری ناشی از وزن 
سازه در نظر گرفته شده است. مدول الاستیسیته 


2.8*10- و جرم متمرکز 1-306.122100 


11 
در نظر گرفته شده اند. جدول (۱) جحگونگی تغییرات 
جرم و سختی سازه را در ارتفاع نشان می دهد. 
به کمک اطلاعات جدول (۱) تابع سختی خمشی 
(*)۳71 < (۶) 1 به صورت زیر درون پابی له اه 
*< 107 2 (10*0 < () و1 


)932 8*2 -1975.9* + 2142.7( 


0 > >1 


(#سه 


(۳( 
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تحلیل پاسخ دینامیکی سازه های بلند تحت تثیر نیروی محوری 


جدول ۱ تغییرات سختی و جرم سازه مثال ۱-۵ در ارتفاع 


1۳ 21-5 1325 235 0 ارتفاع 
تک 7 408.776 | 414388 | 350.102 | جرم واحد طول 
که ۳ ۱۱۵ 11962 7-2-1 2165 ممان اینرسی 
1۳ 16 2245 15+ 2-5 ارتفاع 
ك 2-1 3 362.653 | 381.171 | جرم واحد طول 
9 7 1205.11 | 1333.32 | 1177.18 ممان اینرسی 
جدول ۲ پنج فرکانس اول سازه مثال ۱-۵ 
5 3 0 5 51 
1109 26 119196 2-6 0915 گت 
1071 221290 110۹45 2-939 ۱914 4 
م 
411 221.8 10948 2-938 ۱8914 رک | 3 
ِ 
204 21-38 10983 2-938 ۱8914 ۳-8 3 
2-993 2113 ۱۱0/۵ 2938 0۱914 ۳-0 0 
2-1 2113 7 39.38 ۱8914 ۳-12 
٩۸۳ - 0 6.7759 | 37.21 ۳ 6 18473 ۳ 5‏ 


با توجه به تغییرات ناچیز جرم واحد طول در 
ارتفای جرم واحد طول تیر جایگزین برابر میانگین 
مقادیر )۱( و 


۳ مب < (۳00 در نظر گرفته شده است. 
10 


ارائه شده در حدول برابر 


همچنین تابع نیروی محوری ناشی از وزن جرم متمرکز 
در انتهای آزاد و جرم توزیع شده در ارتفاع سازه به 
صورت زیر در نظر گرفته ننن‌ه. است: 


+3001.83 2 (-1) ,]و + ع <<( 
۸0۱ (-283301.61580 

(۳۵( 
پنج فرکانس اول سازه مطابق روش ارایه شده در 
این مقاله محاسبه شده است. سازه مورد نظر در نرم 
افزار 8۸۳-2000 نیز مدل سازی شده و نتایج تحلیل و 
مدل سازی در نرم افزار در جدول (۲) ارائه شده است. 


جهت در نظر گرفتن اثرات نیروی محوری بر فرکانس 


نشریة مهندسی عمران فردوسی 


طبیعی ارتعاش سازه در تحلیل نرم افزار ابتدا یک آنالیز 
استاتیکی غیر خطی با در نظر گرفتن آثار ۸- ۳انجام 
شده و از ماتریس سختی اصلاح شده در انتهای آتالیز 
۳-۸ جهت انجام آنالیز مودال و محاسبه فرکانس ها 
استفاده شده است. 

طبق رابطه (ع۲)» تابع مد شکل ارتعاش به کمک 
سری توانی ۷ < - (۷ تقریب زده شده است. در 
این روش انتخاب ۳ مناسب به عهده کاربر می باشد. با 
در نظر گرفته می شود که این مساله موجب افزایش 
دقت نتایج تحلیل می شود. به طوری که با در نظر 
گرفتن ۳-3 (در نظر گرفتن چهار جمله اول سری) 
فرکانس های اول و دوم سازه همگرا شده اند که نشان 
از همگرایی مرج روش ارائه شده می باشد. در ازای 
۳-3 فرکانس مد سوم ارتعاش 1796 با جواب تحلیل 
اجزا محدود اختلاف دارد. در حالی که با افزایش ۳ و 


سال سی و یکم شمارة دوه ۱۳۹۷ 


مهرداد محمدنژادسحسن حاج ی کاظمی 


کن نظر گر فین ۳-5 فرکانس مد سوم همگرا شده و این 
اختلاف به 1026 کاهش پیدا کرده است. . فر کانس مد 
چهارم با ۳-8 و فرکانس مد پنجم ارتعاش با ۳-10 
همگرا شده اند. 


سازه بلند با سختی برشی. در این قسمت یک سازه 
بلند ۲۵ طبقه با سیستم قاب محیطی تو در تو که پیش 
از این توسط سایر محققین بررسی شده است به کمک 
روش ارائه شده در این تحقیق تحلیل می گردد. 
ی ی ات 
ترئیسب برایر ]1۵ 3528723010 (ویگ) و 
ت10 7.5538<101 2 ر(و1)در نظر گرفته شده 
است. سختی برشی . جرم واحد طول سازه و ارتفاع 
دق و تسس انس الط 10و ری 
25.8 0 و 0 .۲-75 در نظر گرفته شده 
1 دو فرکانس طبیعی اول ارتعاش سازه محاسبه 
شده و در جدول (۳) با نتایم سایر محققین مقایسه شده 


استت. 


برسی تأثیر تغییر سختی برشی و سختی خمشی بر 
فرکانس طبیعی ارتعاش. در این مثال یک سازه بلند 
۷۰ طبقه با سیستم باربر جانبی قاب محیطی مرکب شده 
با دیوار برشی به عنوان سازه مبنا معرفی می گردد. 
مشخصات سازه ای سیستم باربر جانبی سازه مبنا به 


شرح زیر می باشد: 


1۹ 


سختی برشی معادل 160 *277.56<10 وک: سختی 
خمشی معادل "ه- 1۵ 2.61710 - وک ارتفاع سازه 
2210 11, جرم واحد طول */ *5.ع 2681408 2 
. سه فرکانس اول سازه مبنا مطابق روش ارائه شده در 


این ٌ تحفیو به شرح زیر محاسبه شده است: 


۱ 
900 9۹0 
120 

2 ور 
90 


(۳۹( 
در رابطه (۳۷) زیرنویس "05" به معنای سازه مبنا می 


خمشی. جرم و ارتفاع سازه بر فرکانس های طبیعی 
1 1۳0 وک وک 
() () م(فک) بج( ک) 


بین ۱ تا ۲ تغییر داده می شود. در این نسبت ها روک 


7[ ,10 ریک مشخصات سازه جدید می باشند. فرکانس 
های سازه جدید به صورت 5۵ ,62 ,,(5 در نظر گرفته 
خحمشی. جرم و ارتفاع سازه بر فرکانس طبیعی ارتعاش 
سازه مبنا به صورت زیر محاسبه می شود: 


(۳۷ / سس کت - 1 
/ 


نتایج تحلیل در اشکال (۳) و (4) ارائه شده است. 


حدول ۳ دو فرکانس طبیعی اول ارتعاش سازه مثال ۲-۵ 


سال سی و یکم شماره دی ۱۳۹۷ 


4 ,0011 ۷ 4 متا ۷ 
(0مطامصط ممتاکوم مه متاک ۵0۵ | (0مطامعط امجصممهام015 100) 
2279 2-157 
1921 ۰ 
7 ,1,66 9 ,۷۷۵۸2 8 ,۵8۳16271,] 
3318 32462 315 
20.63 2120 2120 


| روش اراه شد 
3.705 2-056 60: /62۵ 6۵ 
101237 16-1326 (560 /۲۵۵) (62 
60 ,۷۷۵۱۵ 19962 ,۷۷/298 
2-۹61 32462 62۵/60 6۵ 
19239 2-5 (560 ۵07 622 
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که 2 کد 1 که 0 که ‌ 3 ک 
16/)8(۶ ۷5/)6(۶ 
(«لف) (ب) 


شکل ۳ تغییرات فرکانس سازه مبنا بر حسب تغییر پارامترهای سازه ای الف: تغییر سختی خمشی. ب: تغییر سختی برشی 


۱0061 -0-۱0062 3-3 


۱۳/)(9۶ 


(الف) (ب) 


گ 2 ک که 
5ظ(۷5/۱5 , عط(۷۵/)۷۵ 


(پ) 
شکل ۶ تغییرات فرکانس سازه مبنا بر حسب تغییر پارامترهای سازه ای. الف) تغییر همزمان سختی برشی و سختی خمشی. ب) تغییر جرم 


سازه. ‌( تخییر ارتفاع سازه 


نشرية مهندسی عمران فردوسی سال سی و یکم» شمارهٌ دی ۱۳۹۷ 


مهرداد محمدنژادسحسن حاج ی کاظمی 


نتایج شکل (۳-الف) نشان می دهد که تغییسرات 
سختی خمشی یک سازه بلند کمترین تأثیر را بر 
فرکانس مد اول و بیشترین تأثیر را بر فرکانس مد سوم 
ارتعاش سازه مبنا دارد و بالعکس. نتایج شکل (۲-ب) 
نشان می دهد که تغییرات سختی برشی یک سازه 
بیشترین تأثیر را بر مد اول ارتعاش و کمترین تأثیر را بر 
مد سوم ارتعاش سازه مبنا دارد. نتایج شکل (4-الف) 
حاکی از آن می باشد که تغییر همزمان سختی برشی و 
سختی خمشی اثر کاملا یکسان بر تغییر هر سه فرکانس 
اول سازه دارد. نتایج شکل (4-ب) نشان می دهد که 
تغییر جرم سازه نیز تأثیر یکسان بر تغییر هر سه مد اول 
ارتعاش سازه داشته است. و در نهایت نتایج شکل (4- 
پ) نشان می دهد که تغییر ارتفاع سازه بیشترین تأثیر را 
بر تغییر فرکانس سوم و کمترین تأثیر را بر فرکانس اول 


ارتعاش سازه دارد. 


نتیجه گیری 
روش‌های تحلیل اجزا محدود بر پایه مدلسازی گسسته 
استوار است و نیاز به حل همزمان هزاران معادله خطی 
به منظور ارائه نتایج عددی و کاربردی تحلیل سازه دارد. 
همچنین مدلسازی یک سازه بلند در نرم افزار تحلیل 
اجزا محدود بدلیل داشتن المان های سازه ای زیاد؛ 
بسیار وقت گیر می باشد. بنابراین روش تحلیل اجزا 
محدود برای مراحل نهایی و تحلیل دقیق یک سازه بلند 
مناسب می باشد. روش ارائه شده در این تحقیق یک 
دید روشن از تحلیل ارتعاش آزاد سازه بلند ارائه می 
"کنات سادگی و دفت بالای روش ارائه شده در این 
تحقیق موجب می شود که در تحلیل تقریبی و مراحل 
مقدماتی طراحی یک سازه بلند ( که نیازی به مدلسازی 
دقیق نرم افزاری وجود ندارد) قابل استفاده باشد. در این 
تحقیق. سازه بلند در یک محیط پیوسته مدلسازی شده 
ایتشا که آنی ال توافت ان تیوتی نم انیت 
همچنین روش ارائه شده تنها نیاز به حل یک معادله 


نشرية مهندسی عمران فردوسی 


[ 


دیفرانسیل داشته و به راحتی قابل کدنویسی در نرم 
افزارهای مختلف می باشد. نتایج تحلیل تقریبی بدست 
آمده به کمک روش ارائه شده می تواند در جهعت 
کاهش خطاهای مدلسازی دقیق سازه در نرم افزار 
تحلیل اجزا محدود. تخمین ابعاد اولیه المان های سازه 
ای و مشخصات مصالح مورد بررسی و ... مفید واقع 
گردد. 

در این تحقیق. معادلات انتگرالی تضعیف شده 
ماکان رهای ای رتاش اوه بلق 
با خصوصیات سختی و جرم متغیر در ارتفاع تحت تأثیر 
نیروی محوری ارائه شده است. از طریق چهار مرتبه 
انتگرال گیری پی در پی. معادله دیفرانسیل حاکم بر 
ارتعاش تیر به فرم تضعیف شده آن تبدیل شده است. 
تابع مد شکل ارتعاش توسط یک سری توانی تقریب 
زده شده و به این طریق معادلات تضعیف شده تبدیل 
به دستگاه معادلات جبری خحطی شده اند. فرکانس 
طبیعی ارتعاش سازه بلند با محاسبه جواب غیر بدیهی 
دستگاه معادلات محاسبه شده است. نتایج تحلیل نشان 
دهنده همگرایی سریع روش ارایه شده می باشد. به 
طوری که با در نظر گرفتن تنها ۶ جمله از سری توانی 
(۲۶). فرکانس های اول و دوم ارتعاش به جواب تحلیل 
اجزا محدود همگرا شده است. با افزایش تعداد جملات 
در نظر گرفته شده برای سری توانی» دقت فرکانس های 
محاسبه شده افزایش يافته و فرکانس های طبیعی 
مدهای ارتعاشی بالاتر نیز محاسبه شده اند. نتایج تحلیل 
نشان می دهد که تغییرات سختی برشی یک سازه 
بیشترین تأثیر را بر مد اول ارتعاش و کمترین تأثیر را بر 
مد سوم ارتعاش سازه مبنا دارد و بالعکس تغییرات 
سختی خمشی یک سازه بلند کمترین تأثیر را بر 
فرکانس مد اول وبیشترین تأثیر را بر فرکانس مد سوم 
از ان اکتا تعازد: 


سال سی و یکم شمارة دوه ۱۳۹۷ 


اف تحلیل پاسخ دینامیکی سازه های بلند تحت تأثیر نیروی محوری 


مت 


,"صمتاها ین حصبامتاجم0) :۱۷۲۵0۵1 ۱۷۲۵1۲ حمطومحصرور۳ وبتاحصه‌ااض م۸" ,.ظ صم‌لان0 0صه .1 162222 
.(2012) ,324-542 .00 و21 ۷۵۱۰ و.ااظ ,5۵66 1011 «وزکهرز 

٩۳۱۵۱۵۲۵۹‏ مصصم۵۱۱ ۲۷۷ ۵۶ ولوروامصه امتمنطد کصمام)۰۳۶ ر.تنا روعبا 240 .وا متطلکا بکا۷ لته 
(2014) ,740-759 .00 و23 ۷۵۱۰ .با .566 1011 «وزوع(] 

۶ ایوگ 0مصاصاصورمت ۵۶ وتورلفمن متصصموونا وی رهم۱21۵ فصع .2 رط02720امحصطه۱۷ ریک ۲۵2022 
و 5۳۸6۵۲۵ 1۳6 و"فممتامم۲ مصتامی-ظ همتونا ۱۷۲۵۵0 مت‌نعاون ها وله ۲۷۷ توعد 4صره وماب ۲ ۲۲۵۲۵0 
(2014) 10.1002/621.1201۰ :۱۱۵1 رعومز0ابظ آونءعم5 0وه 1011 ]0 

۶۴ ارحص ۱۷۱۵2/109 ۳۲۵۵ 107 ۷۵0۵1 لهم)هصمم۱]2)۵ مامرصصزه خ ,۷۲۸ رتتا0م92202 0صه بک و2827 
۵ ما0 4ص2 وفناع 1 ااوظ رععمن تحمطه مومت 1 0مصصهز۲ ۵۶ عصااوتعدم صماوری 0مصاممصمن 
,4918-30 .00 و36 ۷۵۱ ص۷00 ۱۷۵۵۱۵۴۳۵۵6۵ 0ع:آمرصض ۵۲ مومت تااممم وان همتا۵مط۵ع0) 
.(2012) 

امسدلا مطا عصتانموومی بو ۱۵00 لاحم مامرصص. ه۸ ویک وت022ع۵20 24 ,۷۲ م۷]2101006[24[ 
۰ ,۱۷0۵061۱۱9 ۱۷۵۱6۱۱۵۵6۵ ۸۵۵۱60 0۲ آمدتیول فص اتباظ 1211 ۶ه عمط ۲۷۵0۵ 4ص ۲۲۵۵0۱2۵۵۵1۵5 
.(2012) ,3419-32 .00 ,36 

۶ ولو رام وفوممموم متصصهوررنا بو م۵ 162 تواقصه صم ویک وته۲۵8202 24 ,۳ م۷]2101006[24[ 
,67-80 00۰ و23 ۷۵۱۰ وکوهاتباظ آمزععم5 ۵80 1011 ]0 «وندعن آه6۵ع50 ۲۳6 ,"فعصللاتباظ ال12 بعااطات ]1 
2014(۰) 


۲۷ ۲۵20 ۲۱۵۱2۵0 ملصنا وممنتمنااه ملظ بمرمافتلن۱۲ ۵۶ ولورلممم لقطم۳]0181 رنهک رقالم2 .. 


2013(۰) ,126-43 .00 و22 ۷۵۱۰ .ابا .566 1۵11 «وزءع(] 
0ص اه رلمتت‌ها۳ هو ممتاممتصماهنا تم ممهممومم افمتارآمممه سعل ۳۸ ویک و460[2مصصصهطه]۱۷ 
۱6۳۳۵۵0 0 ۱/۵۵۵6 00 ودا2 6و۱ هناگ رعصصفعظ ۵۴ ومزم‌ممناوم:۳ لمتنتداط لقط01810 1[ 

.(2015) ,655-674 .00 و55 ,۷۵۱ ,لاه 

۸۵0۲020 262ص مه و وتتاممتمطموظ مه .۲ رتتهکلمه م۳ م40[24مصصمصهطه]۱۷ 
39 ۷۵۱۰ ,۳۱9۱۱۵6۲۱۵ ۵0۳0 56۱6۴6۵ ۵۲ ما0 ما۸۵ و عوتاتهعو۳۵ مامامتنه ۷ 100 وصفه‌ظ ۵۴ وتو راهم 
(2014) ,2561-2572 00۰ 

۵ 1۵ 6ط10هنا0ع اجتعماصا مک لمع ۵۶ ممتاهلامم2 مطا ۵ و رلدزمصمصصصه‌طم 4صه ۲ رتتعلهه 
-977 .00 ,16 ,۷۵۱ ,رهظ وا ععوص از ]۵ سول ماک رقمعنطمنطاه لها ۵۶ علور‌امصه صمته۷10 
(2015) ,999 


1076 ]0 20217515 ط۷۱0210 ۳۳۲۲۵۵ ربا رتتاممتمطعوظ مه ۷۲۰ ر2زمملمصمصهطم۱ رب رتتولله . 


وال ۲۱۵۵6۲۱۵۵۳۵1 صنظ کمنممع ۷( 00 ودا۵ع و۱ هانگ رععع۲01 اهتعنه ماطهتنه۷ تعصا قصصنحه‌ه 
2012(۰) ,561-82 .00 ,43 ,۲۷۵۱ 
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